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ABSTRACT 

A comprehensive study of the effect of an unstable tau neutrino on primordial 
nucleosynthesis is presented. The standard code for nucleosynthesis is mod- 
ified to allow for a massive decaying tau neutrino whose daughter products 
include neutrinos, photons, pairs, and/or noninteracting (sterile) daugh- 
ter products. Tau-neutrino decays influence primordial nucleosynthesis in 
three distinct ways: (i) the energy density of the decaying tau neutrino and 
its daughter products affect the expansion rate tending to increase ^He, D, 
and ^He production; (ii) electromagnetic (EM) decay products heat the EM 
plasma and dilute the baryon-to-photon ratio tending to decrease "^He pro- 
duction and increase D and ^He production; and (iii) electron neutrinos and 
antineutrinos produced by tau-neutrino decays increase the weak rates that 
govern the neutron-to-proton ratio, leading to decreased ^He production for 
short lifetimes (< 30 sec) and masses less than about 10 MeV and increased 
^He production for long lifetimes or large masses. The precise effect of a de- 
caying tau neutrino on the yields of primordial nucleosynthesis and the mass- 
lifetime limits that follow depend crucially upon decay mode. We identify 



four generic decay modes that serve to bracket the wider range of possibih- 
ties: tau neutrino decays to (1) sterile daughter products (e.g., I'r ^ + 4>', 
is a very weakly interacting scalar particle); (2) sterile daughter product(s) 
+ daughter products(s) that interacts electromagnetically (e.g., 1/^- z/^+7); 

(3) electron neutrino + sterile daughter product(s) (e.g., z/^ — * z/g + 0); and 

(4) electron neutrino + daughter product (s) that interact electromagneti- 
cally — >■ z/g + e^). Mass-lifetime limits are derived for the four generic 
decay modes assuming that the abundance of the massive tau neutrino is 
determined by its electroweak annihilations. In general, nucleosynthesis ex- 
cludes a tau-neutrino of mass 0.4 MeV — 30 MeV for lifetimes greater than 
about 300 sec. These nucleosynthesis bounds are timely since the current 
laboratory upper bounds to the tau-neutrino mass are around 30 MeV, and 
together the two bounds very nearly exclude a long-lived tau neutrino more 
massive than about 0.4 MeV. Further, our nucleosynthesis bounds together 
with other astrophysical and laboratory bounds exclude a tau neutrino of 
mass 0.4 MeV — 30 MeV of any lifetime that decays with EM daughter prod- 
uct (s). We use our results to constrain the mass times relic abundance of 
a hypothetical, unstable species with similar decay modes. Finally, we note 
that a tau neutrino of mass 1 MeV to 10 MeV and hfetime 0.1 sec — 10 sec 
whose decay products include an electron neutrino can reduce the ^He yield 
to less than that for two massless neutrino species. This fact could be rel- 
evant if the primordial mass fraction of "^He is found to be less than about 
0.23 and can also lead to a modification of the nucleosynthesis bound to the 
number of hght (<^ 1 MeV) neutrino (and other) particle species. 



1 Introduction 



Primordial nucleosynthesis is one of the cornerstones of the hot big-bang 
cosmology. The agreement between the predictions for the abundances of 
D, ^He, '^He and ^Li and their inferred primordial abundances provides the 
big-bang cosmology's earliest, and perhaps most, stringent test. Further, 
big-bang nucleosynthesis has been used to provide the best determination 
of the baryon density || and to provide crucial tests of particle-physics 
theories, e.g., the stringent bound to the number of light neutrino species 

Over the years various aspects of the effect of a decaying tau neutrino 
on primordial nucleosynthesis have been considered [^, |^, |], |], [1^, |TT|, |1^] . 
Each previous study focused on a specific decay mode and incorporated dif- 
ferent microphysics. To be sure, no one study was complete or exhaustive. 
Our purpose here is to consider all the effects of a decaying tau neutrino on 
nucleosynthesis in an comprehensive and coherent manner. In particular, for 
the first time interactions of decay-produced electron neutrinos and antineu- 
trinos, which can be important for lifetimes shorter than 100 sec or so, are 
taken into account. 

The nucleosynthesis limits to the mass of an unstable tau neutrino are cur- 
rently of great interest as the best laboratory upper mass limits [|1^, 31 MeV 
by the ARGUS Collaboration and 32.6 MeV by the CLEO Collaboration^] 
are tantalizingly close to the mass range excluded by nucleosynthesis, approx- 
imately 0.4 MeV to 30 MeV for lifetimes greater than about 300 sec. If the 
upper range of the cosmologically excluded band can be convincingly shown 
to be greater than the upper bound to the mass from laboratory experiments, 
the two bounds together imply that a long-lived tau-neutrino must be less 
massive than about 0.4 MeV. This was the major motivation for our study. 

The effects of a massive, decaying tau neutrino on primordial nucleosyn- 
thesis fall into three broad categories: (i) the energy density of the tau neu- 
trino and its daughter product (s) increase the expansion rate, tending to 
increase ^He, D, and ^He production; (ii) the electromagnetic (EM) plasma 
is heated by the daughter product (s) that interact electromagnetically (pho- 
tons and pairs), diluting the baryon-to-photon ratio and decreasing "^He 

^Both are 95% C.L. mass limits based upon end-point analyses of tau decays to final 
states containing five pions. The CLEO data set contains 113 such decays and the ARGUS 
data set contains 20 such decays |13[ . 
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production and increasing D and '^He production; and (iii) electron neutrino 
(and antineutrino) daughters increase the weak interaction rates that govern 
the neutron-to-proton ratio, leading to decreased ^He production for short 
lifetimes (< 30 sec) and masses less than about 10 MeV and increased ^He 
production for long lifetimes. Decays that take place long after nucleosyn- 
thesis {ti, ~ 10^ sec — 10^ sec) can lead to the destruction of the light elements 



through fission reactions and additional constraints |]T4|, neither of which are 
considered here. 

In terms of the effects on primordial nucleosynthesis there are, broadly 
speaking, four generic decay modes: 

1. Tau neutrino decays to daughter products that are all sterile, e.g., Ur —>■ 
+ (0 is a very weakly interacting boson). Here, only effect (i) comes 
into play. Aspects of this case were treated in Refs. P, ^ |10|, |Tl|, |12 



the very recent work in Ref. |I2| is the most complete study of this 
mode. 

2. Tau neutrino decays to a sterile daughter product (s) plus a daughter 
product(s) that interacts electromagnetically, e.g., u^- — >• z/^ + 7. Here, 
effects (i) and (ii) come into play. This case was treated in Ref. [|], 
though not specifically for a decaying tau neutrino. 

3. Tau neutrino decays into an electron neutrino and sterile daughter 
product (s), e.g., + 4>. Here, effects (i) and (iii) come into play. 
This case was treated in Ref. ||^; however, the interactions of electron 
neutrinos and antineutrinos with the ambient thermal plasma were not 
taken into account. They can be important: The interaction rate of a 
high-energy electron neutrino produced by the decay of a massive tau 
neutrino relative to the expansion rate T/H ~ (mjy/ MeV)( sec/t). 

4. Tau neutrino decays into an electron neutrino and daughter product (s) 
that interact electromagnetically, e.g., Ur — »• + e^. Here, all three 
effects come into play. Aspects of this case were treated in Ref. [|], 
though interactions of electron neutrinos and antineutrinos with the 
ambient thermal plasma were neglected and the i^e-spectrum was taken 
to be a delta function. 

As we shall emphasize more than once, the effect of a tau neutrino of a 
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given mass and lifetime — and therefore limits to its mass/lifetime — depends 
very much upon decay mode. 

While these four generic decay modes serve to bracket the possibihties, 
the situation is actually somewhat more complicated. Muon neutrinos are 
not completely sterile, as they are strongly coupled to the electromagnetic 
plasma down to temperatures of order a few MeV (times of order a fraction 
of a second), and thus can transfer energy to the electromagnetic plasma. 
However, for lifetimes longer than a few seconds, their interactions with the 
electromagnetic plasma are not very significant (see Ref. |jl5[); and so to a 



reasonable approximation muon-neutrino daughter products can be consid- 
ered sterile. Precisely how much electromagnetic entropy is produced and the 
effect of high-energy neutrinos on the proton-neutron interconversion rates 
depend upon the energy distribution of the daughter products and their in- 
teractions with the ambient plasma (photons, pairs, and neutrinos), which 
in turn depends upon the number of daughter products and the decay matrix 
element. 

Without going to extremes, one can easily identify more than ten possible 
decay modes. However, we believe the four generic decay modes serve well 
to illustrate how the nucleosynthesis mass-lifetime limits depend upon the 
decay mode and provide reasonable estimates thereof. In that regard, input 
assumptions, e.g., the acceptable range for the primordial abundances and 
the relic neutrino abundance^ probably lead to comparable, if not greater, 
uncertainties in the precise limits. 

Finally, a brief summary of our treatment of the microphysics: (1) The 
relic abundance of the tau neutrino is determined by standard electroweak 
annihilations and is assumed to be frozen out at its asymptotic value during 
the epoch of nucleosynthesis, thereafter decreasing due to decays only. Be- 
cause we assume that the relic abundance of the tau neutrino has frozen out 
we cannot accurately treat the case of short lifetimes, < (m^/ MeV)~^sec, 
where inverse decays can significantly affect the tau-neutrino abundance and 
that of its daughter products |16[-Q (2) Sterile daughter products, other 
than neutrinos, are assumed to have a negligible pre-decay abundance (if 

^The variation between different calculations of the tau-neutrino abundance are of the 
order of 10% to 20%; they arise to different treatments of thermal averaging, particle 
statistics, and so on. Since we use the asymptotic value of the tau-neutrino abundance 
our abundances are in general smaller, making our limits more conservative. 

•^For generic decay mode (1) the effect of inverse decays for short lifetimes was consid- 
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this is not true, the nucleosynthesis hmits become even more stringent). (3) 
The electromagnetic energy produced by tau-neutrino decays is assumed to 
be quickly thermalized and to increase the entropy in the electromagnetic 
plasma according to the first law of thermodynamics. (4) The perturba- 
tions to the phase-space distributions of electron and muon neutrinos due to 
tau-neutrino decays and partial coupling to the electromagnetic plasma are 
computed. (5) The change in the weak rates that interconvert neutrons and 
protons due to the distorted electron-neutrino distribution are calculated. (6) 
The total energy of the Universe includes that of photons, pairs, neutrinos, 
and sterile daughter product (s). 

The paper is organized as follows; in the next Section we discuss the 
modifications that we have made to the standard nucleosynthesis code. In 
Section 3 we present our results, discussing how a decaying tau neutrino 
affects the yields of nucleosynthesis and deriving the mass/lifetime limits for 
the four generic decay modes. In Section 4 we discuss other astrophysical 
and laboratory limits to the mass/lifetime of the tau neutrino, and finish in 
Section 5 with a brief summary and concluding remarks. 

2 Modifications to the Standard Code 

In the standard treatment of nucleosynthesis it is assumed that there 
are three massless neutrino species that completely decouple from the elec- 
tromagnetic plasma at a temperature well above that of the electron mass 
(T ~ lOMeV ^ me). Thus the neutrino species do not interact with the 
electromagnetic plasma and do not share in the "heating" of the photons 
when the pairs disappear. 

In order to treat the most general case of a decaying tau neutrino we have 
made a number of modifications to the standard code. These modifications 
are of four kinds: (1) Change the total energy density to account for the 
massive tau neutrino and its daughter products; (2) Change the first-law of 
thermodynamics for the electromagnetic plasma to account for the injection 
of energy by tau decays and interactions with the other two neutrino seas; (3) 
Follow the Boltzmann equations for the phase-space distributions for elec- 
tron and muon neutrinos, accounting for their interactions with one another 
and the electromagnetic plasma; (4) Modify the weak interaction rates that 

ered in Ref. it leads to additional mass constraints for short lifetimes. 
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interconvert neutrons and protons to take account of the perturbations to 
the electron-neutrino spectrum. 

These modifications required tracking five quantities as a function of T = 
R~^, the neutrino temperature in the fully decoupled limit {R = the cosmic- 
scale factor). They are: p^^, (where (j) is any sterile, relativistic decay 
product), Ty, and Ag and A^, the perturbations to the electron-neutrino and 
mu-neutrino phase-space distributions. 

Our calculations were done with two separate codes. The first code tracks 
Pu^, p^, Ty, Ag, and A^ as a function of T, for simplicity, using Boltzmann 
statistics. These five quantities were then converted to functions of the pho- 
ton temperature using the T(T^) relationship calculated, and their values 
were then passed to the second code, a modified version of the standard 
nucleosynthesis code ^^T^ .Q We now discuss in more detail the four modifica- 
tions. 



2.1 Energy density 



There are four contributions to the energy density: massive tau neutrinos, 
sterile decay products, two massless neutrino species, and the EM plasma. 
Let us consider each in turn. 

As mentioned earlier, we fix the relic abundance of tau neutrinos assuming 
that freeze out occurs before nucleosynthesis commences (t -C Isec). We 
follow Ref. 



rU| in writing 



(3.151T)2 + m, 
3.151T 



p^{m^ = 0) exp(-t/r,,). 



(1) 



where r is the ratio of the number density of massive neutrinos to a massless 
neutrino species, the (3.151T)^ term takes account of the kinetic energy of 
the neutrino, and the exponential factor takes account of decays. The relic 
abundance is taken from Ref. []10[; for a Dirac neutrino it is assumed that all 
four degrees are freedom are populated for masses greater than 0.3 MeV (see 
Ref. |TD[ for further discussion). 



^The correct statistics for all species are of course used in the nucleosynthesis code; the 
five quantities are passed as fractional changes (to the energy density, temperature and 
rates) to minimize the error made by using Boltzmann statistics in the first code. 
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Note that for temperatures much less than the mass of the tau neutrino, 
Pur/Puimu = 0) = rm,ye~*/'^''/3.151T, which increases as the scale factor 
until the tau neutrinos decay; further, rm^ determines the energy density 
contributed by massive tau neutrinos and hence essentially all of their effects 
on nucleosynthesis. The relic neutrino abundance times mass {ririy) is shown 
in Fig. 1. 

The energy density of the sterile decay products is slightly more compli- 
cated. Since the 0's are massless, their energy density is governed by 

dT T T Ht^' U 
where the first term accounts for the effect of the expansion of the Universe 
and the second accounts for the energy dumped into the sterile sector by 
tau-neutrino decays. The quantity is the fraction of the tau-neutrino 
decay energy that goes into sterile daughters: for p^. all-sterile daughter 
products, ftj) = 1; for ^ (p + I'e or EM, = 0.5; and for all other 
modes = 0. Eq. (||) was integrated numerically, and was passed to the 
nucleosynthesis code by means of a look-up table. 

The neutrino seas were the most complicated to treat. The contribution 
of the neutrino seas was divided into two parts, the standard, unperturbed 
thermal contribution and the perturbation due to the slight coupling of neu- 
trinos to the EM plasma and tau-neutrino decays, 

Pu = PuO + Spu- (3) 

The thermal contribution is simply 6T^/ vr^ per massless neutrino species 
(two in our case). The second term is given as an integral over the pertur- 
bation to the neutrino phase-space distribution, 

'^p-= E Jpd^pMp,t), (4) 

where the factor of two accounts for neutrinos and antineutrinos. 

Finally, there is the energy density of the EM plasma. Since the electro- 
magnetic plasma is in thermal equilibrium it only depends upon T^: 



PEM — 5 1 n 



3K2{me/T- 



1) 



(5) 



where Ki and K2 are modified Bessel functions. We numerically compute 
as a function of T by using the first law of thermodynamics. 
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Pem = -^2' ~^ K2{me/T^). (7) 



2.2 First law of thermodynamics 

Energy conservation in the expanding Universe is governed by the first law 
of thermodynamics, 

C/[PT0T^'] = -pTOTdR\ (6) 

where in our case pxoT = Pem + Pvo + ^Pu + P<t> + Pur, Ptot = Pem + Puo + 
Spu +P4> + Pu-r, SPu = Spu/3, P<j> = P<j>/^, and 

— 1 ^ f 

vr^ vr 

Eq. (P) can be rewritten in a more useful form, 

dXy_ _ -3ij" (ptot + Ptot - 4p^o/3) - d{Sp^ + P<f> + Pur)/dt 
dt dpEu/dTj 

The quantity dpEu/dT^ is easily calculated, and the time derivatives of the 
densities can either be solved for analytically, or taken from the previous 
time step. 

2.3 Neutrino phase-space distribution functions 

The Boltzmann equations governing the neutrino phase-space distribution 
functions in the standard case were derived and solved in Ref. []TB[. We 
briefly summarize that treatment here, focusing on the modifications required 
to include massive tau-neutrino decays. 

We start with the Boltzmann equation for the phase-space distribution 
of neutrino species a in the absence of decays: 



J dU,dU2dU,{27r)^6\pa +P1-P2- Ps) 



dt Ea dE 2Ea processes 

X |A^.+1^2+3r[/a/l-/2/3], (9) 

where the processes summed over include all the standard electroweak 2^2 
interactions of neutrinos with themselves and the electromagnetic plasma, 
and Boltzmann statistics have been used throughout. 

We write the distribution functions for the electron and muon neutrinos 
as an unperturbed part plus a small perturbation: 

/,(p, t) = eM-p/T) + \{p, t), (10) 
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where we have assumed that both species are effectively massless. During 
nucleosynthesis the photon temperature begins to deviate from the neutrino 
temperature T, and we define 

6{t) = TJT - 1. 

Eq. (P) is expanded to lowest order in Aj and leading to master 
equations of the form: 

|;A,(p,t) = 4G|T3[-A,(p,t)A,(p,t) + i?,(p,t)5(t) + Q(p,t) + C:(p,t)], (11) 



where i = e, fi and the expressions for Ai, Bi, Ci, and are given in Ref. |]T5 
[in Eq. (2. lid) for the coefficient (c + 8) should be (c + 7)]. 

In context of tau-neutrino decays we treat decay-produced muon neutri- 
nos as a sterile species, and thus we are only interested in modifying the 
master equation for electron neutrinos to allow for decays. In the case of 
two-body decays (e.g., u^- ^ + (p or Ut- ^ z/e+ EM) the additional term 
that arises on the right-hand side of Eq. ([TT|) is 

|a,(p, T) = --- + ^^S{p- m./2), (12) 

where n^^ is the number density of massive tau neutrinos.0 

The decay mode i/^ — > z/g -|- has a three-body final state, so that the 
energy distribution of electron neutrinos is no longer a delta function. In this 
case, the source term is 

P A f rp-s , 327rV^p(3-4p/mt.) 

-AeP,T =■■■ + — ^p-m^/2), (13 

where for simplicity we have assumed that all particles except the massive 
tau neutrino are ultrarelativistic. 



2.4 Weak- interaction rates 

Given Ae, it is simple to calculate the perturbations to the weak interac- 
tion rates that convert protons to neutrons and vice versa (see Ref. for 

^For TOj^ we actually use our expression for the total tau-neutrino energy = 
^ rrv^ + (3.151T)2. Except for very short lifetimes and small masses, Ey « m^. 
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details). The perturbations to the weak rates are obtained by substituting 
exp(— p/T) + Ae(p, t) for the electron phase-space distribution in the usual 
expressions for the rates |TB[ and then expanding to lowest order. The pertur- 



bations to the rates for proton-to-neutron conversion and neutron-to-proton 
conversion (per nucleon) are respectively 



6Xnp 



-- / EdE{E^ - miy/^E + QfA,{E + Q), (14) 
-- / EdEiE' - mlf/\E - QfK{E - Q), (15) 

AoT„ J me 



where Boltzmann statistics have been used for all species, r„ is the neutron 
mean lifetime, Q = 1.293 MeV is the neutron-proton mass difference, and 



A, 







("^ EdE{E^-mlY'\E-Q)\ 

J me 



The perturbations to the weak rates are computed in the first code and 
passed to the nucleosynthesis code by means of a look-up table. The unper- 
turbed part of the weak rates are computed by numerical integration in the 
nucleosynthesis code; for all calculations we took the neutron mean lifetime 
to be 889 sec. 



3 Results 

In this section we present our results for the four generic decay modes. Mode 
by mode we discuss how the light-element abundances depend upon the mass 
and lifetime of the tau neutrino and derive mass/lifetime limits. We exclude 
a mass and lifetime if, for no value of the baryon-to-photon ratio, the light- 
element abundances can satisfy: 

Yp < 0.24; (16) 
D/H > 10"^; (17) 
(D+='He)/H < 10"^; (18) 
Li/H < 1.4 X 10-^°. (19) 

For further discussion of this choice of constraints to the light-element abun- 
dances we refer the reader to Ref. HI. 
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The ^He and D + ^He abundances play the most important role in deter- 
mining the excluded regions. The mass/lifetime limits that follow necessarily 
depend upon the range of acceptable primordial abundances that one adopts, 
a fact that should be kept in mind when comparing the work of different au- 
thors and assessing confidence levels. Further, the relic abundances used by 
different authors differ by 10% to 20%. Lastly, the precise limit for a specific 
decay mode will of course differ slightly from that derived for its "generic 
class." 

In illustrating how the effects of a decaying tau neutrino depend upon 
lifetime and in comparing different decay modes we use as a standard case 
an initial (i.e., before decay and annihilations) baryon-to-photon ratio 
Tji = 8.25 X 10~^°. In the absence of entropy production (no decaying tau 
neutrino or decay modes 1 and 3 which produce no EM entropy) the final 
baryon-to-photon ratio r^o — 477^/11 = 3x10"^'', where 4/11 is the usual factor 
that arises due to the entropy transfer from pairs to photons. In the case 
of decay modes 2 and 4 there can be significant EM entropy production, and 
the final baryon-to-photon ratio i] = Tjo/^Sf/Si) < rjo {Sf/Si is the ratio of 
the EM entropy per comoving volume after decays to that before decays). 
Even though rjo does not correspond to the present baryon-to-photon ratio 
if there has been entropy production, we believe that comparisons for fixed 
r]o are best for isolating the three different effects of a decaying tau neutrino 
on nucleosynthesis. For reference, in the absence of a decaying tau neutrino 
the ^He mass fraction for our standard case is: Yp — 0.2228 (two massless 
neutrino species) and 0.2371 (three massless neutrino species). 

3.1 z/t- — ^ sterile daughter products 

Since we are considering lifetimes greater than 0.1 sec, by which time muon 
neutrinos are essentially decoupled, the muon neutrino is by our definition 

effectively sterile, and examples of this decay mode include, z/,- — z/^ + 
where is some very weakly interacting scalar particle (e.g., majoron) or 

l^T ^ i^n + i^n + ^n- 

For this decay mode the only effect of the unstable tau neutrino on nucle- 
osynthesis involves the energy density it and its daughter products contribute. 
Thus, it is the simplest case, and we use it as "benchmark" for comparison 
to the other decay modes. The light-element abundances as a function of 
tau-neutrino lifetime are shown in Figs. 2-4 for a Dirac neutrino of mass 
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20 MeV. 

The energy density of the massive tau neutrino grows relative to a mass- 
less neutrino species as rm^jZT until the tau neutrino decays, after which 
the ratio of energy density in the daughter products to a massless neutrino 
species remains constant. For tau-neutrino masses in the 0.3 MeV to 30 MeV 
mass range and lifetimes greater than about a second the energy density of 
the massive tau neutrino exceeds that of a massless neutrino species before 
it decays, in spite of its smaller abundance (i.e., r ^ 1). The higher energy 
density increases the expansion rate and ultimately ^He production because 
it causes the neutron-to-proton ratio to freeze out earlier and at a higher 
value and because fewer neutrons decay before nucleosynthesis begins. Since 
the neutron-to-proton ratio freezes out around 1 sec and nucleosynthesis oc- 
curs at around a few hundred seconds, the ^He abundance is only sensitive 
to the expansion rate between one and a few hundred seconds. 

In Fig. 2 we see that for short lifetimes [jy ^ 1 sec) the '^Hc mass fraction 
approaches that for two massless neutrinos (tau neutrinos decay before their 
energy density becomes significant). As expected, the "^He mass fraction 
increases with lifetime leveling off at a few hundred seconds at a value that 
is significantly greater than that for three massless neutrino species. 

The yields of D and ^Hc depend upon how much of these isotopes are not 
burnt to "^He. This in turn depends upon competition between the expansion 
rate and nuclear reaction rates: Faster expansion results in more unburnt D 
and ^He. Thus the yields of D and ^He increase with tau-neutrino hfetime, 
and begin to level off for lifetimes of a few hundred seconds as this is when 
nucleosynthesis is taking place (see Fig. 3). 

The effect on the yield of ^Li is a bit more complicated. Lithium produc- 
tion decreases with increasing r/ for 77 < 3 x 10~^° because the final abundance 
is determined by competition between the expansion rate and nuclear pro- 
cesses that destroy ^Li, and increases with increasing 77 for > 3 x 10~^° 
because the final abundance is determined by competition between the ex- 
pansion rate and nuclear processes that produce ''Li. Thus, an increase in 
expansion rate leads to increased ''Li production for 77 < 3 x 10^^° and de- 
creased ^Li production for 77 > 3 x 10~^°; this is shown in Fig. 4. Put another 
way the valley in the ^Li production curve shifts to larger r] with increasing 
tau-neutrino lifetime. 

We show in Figs. 5 and 6 the excluded region of the mass/lifetime plane 
for a Dirac and Majorana tau neutrino respectively. As expected, the ex- 
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eluded mass range grows with lifetime, asymptotically approaching 0.3 MeV 
to 33 MeV (Dirac) and 0.4 MeV to 30 MeV (Majorana). We note the signif- 
icant dependence of the excluded region on hfetime; our results are in good 
agreement with the one other work where comparison is straightforward [|l^ 
and in general agreement with Refs. R, |12 



3.2 Vj- sterile + electromagnetic daughter products 

Again, based upon our definition of sterility, the sterile daughter could be a 
muon neutrino; thus, examples of this generic decay mode include z/^ — > ^^^+7 
or Ur ^ + e^. Our results here are based upon a two-body decay (e.g., 
i^T — * + l), and change only slightly in the case of a three-body decay 
(e.g., Vt ^ T^ii + e^), where a larger fraction of the tau-neutrino mass goes 
into electromagnetic entropy. 

Two effects now come into play: the energy density of the massive tau 
neutrino and its daughter products speed up the expansion rate, tending 
to increase ^He, ^He, and D production; and EM entropy production due 
to tau-neutrino decays reduce the baryon-to-photon ratio (at the time of 
nucleosynthesis), tending to decrease "^He production and to increase D and 
^He production. Both effects tend to shift the ^Li valley (as a function of ?7o) 
to larger r]Q. 

While the two effects have the "same sign" for D, ^He, and "^Li, they 
have opposite signs for "^He. It is instructive to compare ^He production as a 
function of hfetime to the previous "all-sterile" decay mode. Because of the 
effect of entropy production, there is little increase in "^He production until 
a lifetime greater than 1000 sec or so. For lifetimes greater than 1000 sec the 
bulk of the entropy release takes place after nucleosynthesis, and therefore 
does not affect the value of i] during nucleosynthesis. 

Because of the competing effects on ^He production, the impact of an 
unstable, massive tau neutrino on nucleosynthesis is significantly less than 
that in the all-sterile decay mode for lifetimes less than about 1000 sec. The 
excluded region of the mass/lifetime plane is shown in Figs. 5 and 6. For 
lifetimes greater than about 1000 sec the excluded mass interval is essentially 
the same as that for the all-sterile decay mode; for shorter lifetimes it is 
significantly smaller. 

Finally, because of entropy production, the final value of the baryon-to- 
photon ratio is smaller for fixed initial baryon-to-photon ratio: it is reduced 
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by the factor by which the entropy per comoving volume is increased. In the 
hmit of significant entropy production {Sf/Si ^ 1), this factor is given by, 
cf. Eq. (5.73) of Ref. |g, 



I TTn I T 

VS.-0,13r,„„^.l,5^;^. (20) 

A precise calculation of entropy production for this decay mode is shown in 
Fig. 7. As can be seen in the figure or from Eq. (PD|), entropy production 
becomes significant for lifetimes longer than about 100 sec. 

3.3 Ur ^ i^e -\- sterile daughter products 

Once again, by our definition of sterility this includes decay modes such as 
z/^ — > i/g + or z/,- — > z/g + z/^z/^. Here, we specifically considered the two-body 
decay mode z/,- — z/g + 0, though the results for the three-body mode are 
very similar. 

Two effects come into play: the energy density of the massive tau neutrino 
and its daughter products and the interaction of daughter electron neutrinos 
with the nucleons and the ambient plasma. The first effect has been discussed 
previously. The second effect leads to some interesting new effects. 

Electron neutrinos and antineutrinos produced by tau-neutrino decays 
increase the weak rates that govern the neutron-to-proton ratio. For short 
lifetimes (< 30 sec) and masses less than about 10 MeV the main effect is 
to delay slightly the "freeze out" of the neutron-to-proton ratio, thereby 
decreasing the neutron fraction at the time of nucleosynthesis and ultimately 
■^He production. For long lifetimes, or short lifetimes and large masses, the 
perturbations to the n — > p and p n rates (per nucleon) are comparable; 
since after freeze out of the neutron-to-proton ratio there are about six times 
as many protons as neutrons, this has the effect of increasing the neutron 
fraction and ^He production. This is illustrated in Fig. 8. The slight shift 
in the neutron fraction does not affect the other light-element abundances 
significantly. 

The excluded portion of the mass/lifetime plane is shown in Figs. 5 and 
6. It agrees qualitatively with the results of Ref. [§.0 Comparing the limits 

^Thc authors of Ref. use a less stringent constraint to ^He production, Yp < 0.26; 
in spite of this, in some regions of the m,y — Ti, plane their bounds are as, or even more, 
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for this decay mode with the all-sterile mode, the effects of electron-neutrino 
daughter products are clear: for long lifetimes much higher mass tau neutri- 
nos are excluded and for short lifetimes low-mass tau neutrinos are allowed. 

3.4 Ur ^ i^e -\- electromagnetic daughter products 

Now we consider the most complex of the decay modes, where none of the 
daughter products is sterile. Specifically, we consider the decay mode — *• 
z/g -|- e^, though our results change very little for the two-body decay Uj. — *■ 

+ 

In this case all three effects previously discussed come into play: the 
energy density of the massive tau neutrino and its daughter products speed 
up the expansion rate; the entropy released dilutes the baryon-to-photon 
ratio; and daughter electron neutrinos increase the weak- interact ion rates 
that control the neutron fraction. The net effect on ^He production is shown 
in Fig. 9 for a variety of tau-neutrino masses. The main difference between 
this decay mode and the previous one, ^ + sterile, is for lifetimes 
between 30 sec and 300 sec, where the increase in "^He production is less due 
to the entropy production which reduces the baryon-to-photon ratio at the 
time of nucleosynthesis. 

The excluded region of the mass/lifetime plane is shown in Figs. 5 and 6. 
It agrees qualitatively with the results of Ref . .[] The excluded region for 
this decay mode is similar to that of the previous decay mode, except that 
lifetimes less than about 100 sec are not excluded as entropy production has 
diminished ^He production in this lifetime interval. 

3.5 Limits to a generic light species 

We can apply the arguments for the four decay modes discussed above to a 
hypothetical species whose relic abundance has frozen out at a value r relative 
to a massless neutrino species before the epoch of primordial nucleosynthesis 

stringent. This is presumably due to the neglect of electron-neutrino interactions with the 
ambient plasma. 

^The authors of Ref. use a less stringent constraint to ^He production, Yp < 0.26; 
in spite of this, in some regions of the rrii, — plane their bounds are as, or even more, 
stringent. This is presumably due to the neglect of electron- neutrino interactions with the 
ambient plasma. 



14 



(also see Refs. [§|, 0]). The previous limits become limits to rm as a function 
of lifetime r and mass m, which are difficult to display. With the exception 
of the effect that involves daughter electron neutrinos, all other effects only 
depend upon rm, which sets the energy density of the massive particle and 
its daughter products. In Fig. 10, we show that for lifetimes greater than 
about 100 sec and masses greater than about 10 MeV, the ^He production is 
relatively insensitive to the mass of the decaying particle. This means that 
for lifetimes greater than about 100 sec the limit to rm should be relatively 
insensitive to particle mass. 

We show in Fig. 11 the excluded regions of the rm-r plane for a 20 MeV 
decaying particle. In deriving these limits we used the same criteria for 
acceptable light-element abundances and assumed three massless neutrino 
species. The limits to rm for decay modes without electron-neutrino daughter 
products are strictly independent of mass; the two other should be relatively 
insensitive to the particle mass for r > 100 sec (and the actual limits are 
more stringent for m > 20 MeV). 



4 Laboratory and Other Limits 

There are a host of other constraints to the mass and lifetime of the tau 
neutrino ||2^. As a general rule, cosmological arguments, such as the one 
presented above, pose upper limits to the tau-neutrino lifetime for a given 
mass: cosmology has nothing to say about a particle that decays very early 
since it would not have affected the "known cosmological history." Labora- 
tory experiments on the other hand pose lower limits to the lifetime because 
nothing happens inside a detector if the lifetime of the decaying particle is 
too long. Finally, astrophysical considerations generally rule out bands of 
lifetime since "signals" can only be detected if (a) the tau neutrinos escape 
the object of interest before decaying and (b) decay before they pass by 
earthly detectors. 

4.1 Laboratory 

The most important limits of course are the direct limits to the tau-neutrino 
mass. These have come down steadily over the past few years. The current 
upper limits are 31 MeV and 32.6 MeV [|13 . 
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If the tau neutrino has a mass greater than 2me = 1.02 MeV, then the 
decay i/,- — > z/g + takes place through ordinary electroweak interactions at 
a rate 

r = 9>i\u \'\U P ^ K/MeV)^|^,.p 

1927r3 ' ' 2.9 X 10^ sec ' ^ ^ 

where t/gr and ?7ee are elements of the unitary matrix that relates mass eigen- 
states to weak eigenstates, the leptonic equivalent of the Cabbibo-Kobayashi- 
Maskawa matrix. We note that the rate could be larger (or even perhaps 
smaller) in models where the decay proceeds through new interactions. Thus, 
limits to Uer give rise to model-dependent limits to the tau-neutrino lifetime. 

A number of experiments have set limits to Uer- The most sensitive 
experiment in the mass range 1.5 MeV < rriu < 4 MeV was performed at the 
power reactor in Gosgen, Switzerland which produces tau neutrinos at a 
rate proportional to |f/erP through decay of heavy nuclei and z/g — Ut mixing. 
Above this mass range, experiments that search for additional peaks in the 
positron spectrum of the vr"*" ^ e"*"// decay (due to — Vr mixing) provide 
the strictest limits. In the mass range 4 MeV < my < 20 MeV, Bryman et 



al. [22| set the limits shown in Fig. 12; for larger masses the best limits come 



from Ref. [|3|. 



There are also direct accelerator bounds to the lifetime of a an unsta- 
ble tau neutrino that produces a photon or pair. In particular, as has 
been recently emphasized by Babu et al. [^], the BEBC beam dump exper- 



iment p5[ provides model-independent limits based upon the direct search 
for the EM decay products. These limits, while not quite as strict as those 
mentioned above, are of interest since they apply to the photon mode and to 
the mode even if the decay proceeds through new interactions. The limit, 

> 0. 18 [mj MeV) sec, (22) 

is shown in Fig. 12. 



4.2 Astrophysical 

The standard picture of type II supernovae has the binding energy of the 
newly born neutron star (about 3 x 10^^ erg) shared equally by neutrinos of 
all species emitted from a neutrinosphere of temperature of about 4 MeV. 
There are two types of limits based upon SN 1987A, and combined they rule 
out a large region of rriy — plane. 
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First, if the tau neutrino decayed after it left the progenitor supergiant, 
which has a radius i? ~ 3 x 10^^ cm, the high-energy daughter photons 
could have been detected pB| , p7| The Solar Maximum Mission (SMM) 



Gamma-ray Spectrometer set an upper limit to the fluence of 7 rays during 
the ten seconds in which neutrinos were detected: 

< 0.9 cm-2; 4.1 MeV < < 6.4 MeV. (23) 

As we will see shortly, if only one in lO^'' of the tau neutrinos leaving the 
supernova produced a photon, this limit would have been saturated. In 
the mass regime of interest there are two ways out of this constraint: The 
lifetime can be so long that the arrival time was more than ten seconds after 
the electron antineutrinos arrived, or the lifetime can be so short that the 
daughter photons were produced inside the progenitor. We can take account 
of both of these possibilities in the following formula for the expected fluence 
of 7 rays: 

A,io = UpW^B^{F,F2) (24) 

where the subscript 10 reminds us that we are only interested in the first ten 
seconds, fi^p ~ 1.4 x 10^° cm~^ is the fluence of a massless neutrino species, 
Wj ~ 1/4 is the fraction of decay photons produced with energies between 
4.1 MeV and 6.4 MeV, Fi is the fraction of tau neutrinos that decay outside 
the progenitor, and F2 is the fraction of these that decay early enough so that 
the decay products were delayed by less than ten seconds. The quantity is 
the branching ratio to a decay mode that includes a photon. For m,^ > 1 MeV 
one expects the z/g -|- mode to be dominant; however, ordinary radiative 
corrections should lead to ~ 10~^ [|^. Finally angular brackets denote 
an average over the Fermi-Dirac distribution of neutrino momenta, 

^ 1.5C(3)r3/o e^/+l' ^^^^ 
where T ~ 4 MeV is the temperature of the neutrinosphere and E = {p^ + 

To evaluate the fluence of gamma rays we need to know Fi and F2. The 
fraction Fi that decay outside the progenitor is simply e"*^/^^ where ti = 
R/v = RE/p and the "lab" lifetime tl = rE/rriy. Of these, the fraction 
whose decay products arrive after ten seconds is e~*2/^^/e"*^/^^ where ^2 = 
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10sec/(l — f/c); thus, F2 = I — e^*!^*^)/^!,^ Figure 12 shows this constraint 
assuming a branching ratio = 10^^. 

The second constraint comes from observing that if tau neutrinos de- 
cayed within the progenitor supergiant, the energy deposited (up to about 
]^g53 gj,g-j -\yould have "heated up" the progenitor so much as to conflict with 
the observed optical luminosity of SN 1987A (and other type II supernovae) 



30|] . We require 

^input = ((1 - F,))E, < 10^^ erg, (26) 

where ~ 10^'^ erg is the energy carried off by a massless neutrino species, 
and 1 — -Fi is the fraction of tau neutrinos that decay within the progenitor. 
This constraint is mode-independent since decay-produced photons or 
pairs will equally well "overheat" the progenitor. As Fig. 12 shows, the 
"supernova-light" bound is extremely powerful. 

Finally, a note regarding our SN 1987A constraints. We have assumed 
that a massive tau-neutrino species has a Fermi-Dirac distribution with the 
same temperature as a massless {m^ -C 10 MeV) neutrino species. This is 
almost certainly false. Massive {m^ > 10 MeV or so) tau neutrinos will drop 
out of chemical equilibrium (maintained by pair creation/annihilations and 
possibly decays/inverse decays) interior to the usual neutrinosphere as the 
Boltzmann factor suppresses annihilation and pair creation rates relative to 
scattering rates. This leads us to believe that we have actually underestimated 
the fluence of massive neutrinos. While the problem has yet to be treated 
rigorously, we are confldent that, if anything, our simplifled treatment results 
in limits that are overly conservative. Accurate limits await a more detailed 
analysis [E^ . 



4.3 Cosmological 

The most stringent cosmological constraint for masses 0.1 MeV < m < 
100 MeV is the nucleosynthesis bound discussed in this paper. Nonetheless, 
it is worthwhile to mention some of the other cosmological limits since they 
are based upon independent arguments. A stable tau neutrino with mass in 
the MeV range contributes much more energy density than is consistent with 
the age of the Universe. Such a neutrino must be unstable, with a lifetime 
short enough for its decay products to lose enough most of their energy to 
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"red shifting " |33]. The hfetime hmit is mass dependent; a neutrino with a 



mass of about 1 MeV must have a hfetime shorter than about 10 sec, and 
the constraint gets less severe for larger or smaller masses. There is an even 
more stringent bound based the necessity of the Universe being matter domi- 
nated by a red shift of about 10^ in order to produce the observed large-scale 
structure |3^. Finally, there are other nucleosynthesis bounds based upon 
the dissociation of the light elements by decay-produced photons or electron- 
neutrinos and by pairs produced by the continuing annihilations of 



tau neutrmos 



5 Summary and Discussion 

We have presented a comprehensive study of the effect of an unstable tau 
neutrino on primordial nucleosynthesis. The effects on the primordial abun- 
dances and the mass/lifetime limits that follow depend crucially upon the 
decay mode. In the context of primordial nucleosynthesis we have identified 
four generic decay modes that bracket the larger range of possibilities: (1) 
all-sterile daughter products; (2) sterile daughter product (s) + EM daughter 
product (s); (3) z/g + sterile daughter product (s); and (4) z/g + EM daughter 
product (s). The excluded regions of the tau-neutrino mass/lifetime plane for 
these four decay modes are shown in Figs. 5 (Dirac) and 6 (Majorana). 

In the limit of long lifetime {t^, ^ 100 sec), the excluded mass range is: 
0.3 MeV - 33 MeV (Dirac) and 0.4 MeV - 30 MeV (Majorana). Together 
with current laboratory upper mass limits, 31 MeV (ARGUS) and 32.6 MeV 
(CLEO), our results very nearly exclude a long-lived, tau neutrino more 
massive than about 0.4 MeV. Moreover, other astrophysical and laboratory 
data exclude a tau-neutrino in the 0.3 MeV — 50 MeV mass range if its decay 
product (s) include a photon or pair. Thus, if the mass of the tau neutrino 
is the range 0.4 MeV to 30 MeV, then its decay products cannot include a 
photon or an pair and its lifetime must be shorter than a few hundred 
seconds. 



We note that the results of Ref. 12 for the all-sterile decay mode are 



more restrictive than ours, excluding masses from about 0.1 MeV to about 
50 MeV for t^, ^ 100 sec. This traces in almost equal parts to (i) small 
{AY ~ +0.003), but significant, corrections to the ^He mass fraction and 
(ii) slightly larger relic neutrino abundance. With regard to the first differ- 
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ence, this illustrates the sensitivity to the third significant figure of the ^He 
mass fraction. With regard to the second difference, it is probably correct 
that within the assumptions made the tau-neutrino abundance during nu- 
cleosynthesis is larger than what we used. However, other effects that have 
been neglected probably lead to differences in the tau-neutrino abundance of 
the same magnitude. For example, for tau-neutrino masses around the upper 
range of excluded masses, 50 MeV — 100 MeV, finite-temperature corrections, 
hadronic final states (e.g., a single pion), and tau-neutrino mixing have not 
been included in the annihilation cross section and are likely to be important 
at the 10% level. 

So is a tau neutrino with lifetime greater than a few hundred seconds and 
mass greater than a fraction of an MeV ruled out or not? Unlike a limit based 
upon a laboratory experiment, it is impossible to place standard error fiags 
on an astrophysical or cosmological bound. This is because of assumptions 
that must be made and modeling that must be done. For example, the 
precise hmits that one derives depend upon the adopted range of acceptable 
light-element abundances. To be specific, in Ref. |jl2|] the upper limit of the 
excluded mass range drops to around 38 MeV and the lower limit increases 
to about 0.4 MeV when the primordial ^He mass fraction is allowed to be 
as large as 0.245 (rather than 0.240). In our opinion, a very strong case 
has been made against a tau-neutrino mass in the mass range 0.4 MeV to 
30 MeV with lifetime much greater than 100 sec; together with the laboratory 
limits this very nearly excludes a long-lived, tau neutrino of mass greater than 
0.4 MeV. 

Perhaps the most interesting thing found in our study is the fact that a tau 
neutrino of mass 1 MeV to 10 MeV and lifetime 0.1 sec to 10 sec that decays 
to an electron neutrino and a sterile daughter product can very significantly 
decrease the ^He mass fraction (to as low as 0.18 or so). It has long been 
realized that the standard picture of nucleosynthesis would be in trouble if 
the primordial ^He mass fraction were found to be smaller than about 0.23; 
within the standard framework we have found one way out: an unstable 
tau neutrino.0 In principle, the possibility of an unstable tau neutrino also 
loosens the primordial-nucleosynthesis bound to the number of light species 

^Based upon dimensional considerations the lifetime for the mode Vr ^ i^e + 4* 
expected to be ^ Sirf'^/mf,, where / is the energy scale of the superweak interactions 
that mediate the decay. For Ti, ~ 10 sec and nii, ~ 10 MeV, / ~ 10^ GeV. 
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which is largely based on the overproduction of ^He. However, an unstable 
tau neutrino does not directly affect the primordial important nucleosynthesis 
bound to the baryon-to- photon ratio (and Q^) as this bound involves the 
abundances of D, ^He, and ^Li and not ^He. 

Finally, we translated our results for the tau neutrino into limits to the 
relic abundance of an unstable, hypothetical particle species that decays into 
one of the four generic decay models discussed. Those very stringent limits 
are shown in Fig. 11. 
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6 Figure Captions 



Figure 1: The relic neutrino abundance used in our calculations as a func- 
tion of neutrino mass: Dirac neutrino, including both helicity states (solid 
curve), and Majorana neutrino (broken curve); results from Ref. [|10|- 

Figure 2: The ^He yield as a function of tau-neutrino lifetime for the four 
generic decay modes, a 20 MeV Dirac neutrino, and a baryon-to-photon ratio 
that in the absence of entropy production leads to a present value rjo = 3 x 
10~^°. For reference, the ^He yield in the absence of a decaying tau neutrino 
is Yp = 0.2228 (2 massless neutrinos) and 0.2371 (3 massless neutrinos). 

Figure 3: The D (solid) and D + ^He (broken) yields as a function of tau- 
neutrino lifetime for the all-sterile decay mode, t] = 3 x lO^^*', and a 20 MeV 
Dirac neutrino. 

Figure 4: The ^Li yield as a function of tau-neutrino lifetime for the all- 
sterile decay mode, 1] = 10~^^ (solid) and 10~^ (broken), and a 20 MeV Dirac 
neutrino. 

Figure 5: Excluded regions of the mass-lifetime for a Dirac neutrino and 
the four generic decay modes. The excluded regions are to the right of the 
curves; our results are not applicable to the region labeled N/A as tau- 
neutrino inverse decays can be important and have not been included (see 
Ref. 0). 

Figure 6: Excluded regions of the mass-lifetime for a Majorana neutrino 
and the four generic decay modes. The excluded regions are to the right of 
the curves; our results are not applicable to the region labeled N/A as inverse 
decays can be important. 

Figure 7: Entropy production as a function of tau-neutrino lifetime for 
a Dirac neutrino of mass 1, 5, 10, 20 MeV and the Ur ^ 4> + EM decay 
mode. Sf/Si is the ratio of the entropy per comoving volume before and 
after tau-neutrino decays. 

Figure 8: ^He yield as a function of tau-neutrino lifetime for the i^r 

Ue + (p decay mode, rj = 3 x 10~^°, and Dirac masses of 1, 5, 10, 20 MeV. 
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For reference, the "^He yield in the absence of a decaying tau neutrino is 
Yp = 0.2228 (2 massless neutrinos) and 0.2371 (3 massless neutrinos). 

Figure 9: ''He yield as a function of tau-neutrino lifetime for the z/^ v^. 
+ EM decay mode, ?7o = 3 x 10~^°, and Dirac masses of 1, 5, 10, 20MeV. 
For reference, the "^He yield in the absence of a decaying tau neutrino is 
Yp = 0.2228 (2 massless neutrinos) and 0.2371 (3 massless neutrinos). 

Figure 10: ^He yield as a function of lifetime for the + </> decay mode, 
r]Q = 3 X 10~^°, rm = 3.5, and masses of 5, 10, 20, 30MeV. For reference, 
the ^He yield in the absence of a decaying tau neutrino is Yp = 0.2228 (2 
massless neutrinos) and 0.2371 (3 massless neutrinos). 

Figure 1 1 : Excluded regions of the rm — t plane for the four different decay 
modes and a 20 MeV mass particle. The limits for the first two decay modes 
are strictly independent of mass; for the last two decay modes they should 
be relatively insensitive to mass for r > 100 sec (and actually more stringent 
than those shown here for m > 20 MeV). Excluded regions are above the 
curves. 

Figure 12: Regions of the tau-neutrino mass-lifetime plane excluded by lab- 
oratory experiments and astrophysical arguments. The excluded regions are 
to the side of the curve on which the label appears. The dashed curve sum- 
marizes a host of different laboratory limits to |t/erP, translated to a model- 
dependent bound to Ti^^Ur — > t'eC^), cf. Eq. (^l]). SNL denotes the supernova- 
light constraint which extends to lifetimes shorter than those shown here. 
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